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I.  INTRODUCTION 


This  report  develops  some  approximation  techniques  for  use  in 
communications  network  design  problems.  The  techniques  deal  with 
"average  properties"  (to  be  defined  later)  of  networks  and  thus  lose 
sight  of  the  discrete  structure;  hence,  they  are  not  suitable  for  such 
problems  as  efficiently  sizing  each  link  in,  say,  the  AUTOVON  network. 
The  techniques  can,  however,  tell  us  sometljing  about  the  average  link 
size,  among  other  things.  Average  rather  than  precise  parameters  are 
already  widely  utilized  in  the  Erlang  formulas  and  in  queueing  theory. 
This  report  is  useful  for  looking  at  average  properties  at  a higher  level 
of  system  issues  in  order  to  provide  a better  feel  for  the  gross  charac- 
teristics of  various  architectural  issues. 

Since  the  techniques  developed  work  with  average  properties  of 
the  network,  they  require  only  that  the  designer  know  average  quantities 
associated  with  the  exogenous  variables  in  network  design.  These  too 
will  be  defined  later.  At  this  point  we  merely  set  the  tone  of  this 
paper;  viz.,  we  wish  to  find  analytical  relationships  between  basic 
structural  properties  (on  the  average)  of  a communications  network  and 
the  external  variables  which  impact  the  network. 

The  development  of  these  techniques  was  motivated  by 
the  need  for  a mechanism  to  assess  quantitatively  various  future  DCS 
architectural  alternatives.  There  are  two  fundamental  ways  of 
approaching  this  problem.  One  is  to  attempt  the  construction  of  network 
design  models  in  the  traditional  sense.  Such  models  require 
detailed  predictions  concerning  the  values  of  external  variables 
at  some  point  in  future  time.  For  example,  one  needs  to  predict  the 
location  and  associated  traffic  volume,  by  type,  of  each  future  DCS 
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user.  This  approach  is  discussed  in  reference  [1]  and  will  not  be 
considered  herein.  The  second  approach  is  the  subject  of  this  paper.  Here 
we  take  a macroscopic  look  at  the  network  and  its  environment.  Only  one 
architectural  issue  is  dealt  with  in  this  paper.  This  is  network  topology  - 
i.e.,  its  effect  on  Network  Grade  of  Service  (GOS)  and  survivability  under 
constrained  costs  in  a terrestrial  backbone  network. 

1.  PROBLEM  DISCUSSION 

A communications  network^ is  to  be  designed  to  service  a rectangular 
area  whose  dimensions  are  A and  B miles,  as  shown  in  Figure  1. 

Very  little  is  known  about  the  requirements  of  the  users.  Consequently, 
using  the  Laplacian  assumption  of  rationality  * [2],  we  assume  that: 

Al.  users  are  uniformly  distributed  in  the  area  and 
all  users  are  equally  likely  to  generate  traffic. 

A2-  A randomly  selected  user  is  equally  likely  to  call  any 
other  user  in  the  area. 

A3.  The  total  offered  traffic  for  the  network  is  E erlangs. 

For  our  purposes  it  is  reasonable  to  assume  that  the  network  is 
constructed  from  two  basic  elements;  nodes  and  links.  A cost  model  for 
each  element  is  given  by  the  assumptions: 

A4.  Link  costs,  D|^,  are  a function  of  the  number  of  channels, 
c,  in  the  link;  and  the  length  in  miles,  I,  of  the  link: 

= kcl 

where  fe  is  the  cost  per  channel  mile.  All  links  are  assumed  to  be 

full  duplex;  where  a "link"  is  taken  to  comprise  the  total  interconnection 

of  two  nodes. 

* This,  in  essence,  means  that  when  we  do  not  know  the  value  of  a para- 
meter we  assume  all  cases  to  occur  with  equal  probability. 


A5.  Node  costs,  Dj^,  comprise  a fixed  cost,  a,  and  a cost 
per  channel  termination,  fa: 

= a + fa-t 

where  t is  the  number  of  channels  terminating  on  the  node. 

We  are  especially  concerned  with  the  design  of  survivable  networks; 
that  is,  networks  which  operate  at  some  minimum  level  after  an  adversary 
has  destroyed  a number  of  nodes.  From  this  it  follows  that  we  restrict 
our  consideration  to  networks  having  a high  degree  of  symmetry;  i.e., 
any  randomly  selected  node  is  connected  to  its  neighbors  in  a uniform 
way  (often  referred  to  as  distributed  networks).  In  this  fashion  the 
enemy  is  given  no  information  as  to  what  part  of  the  network  is  most 
sensitive  - in  fact,  ideally  no  one  part  is  more  sensitive  than  any 
other.  As  a consequence: 

A6.  When  the  network  is  attacked,  nodes  are  removed  by  a 
random  process. 

Our  approach  to  this  problem  requires  that  we  impose  some  geometric 
regularity  on  potential  network  designs  while  leaving  as  design 
variables  the  number  of  nodes  and  the  number  of  links.  The  necessary 
geometric  regularity  consists  of  postulating  a square  grid  arrangement 
of  nodes  in  the  area.  Thus,  if  we  were  to  investigate  a design  con- 
sisting of  M'N  nodes  they  would  be  arranged  in  a uniform  grid  inside 
the  area, as  shown  in  Figure  2.  We  assume,  for  the  sake  of  definiteness: 

A7.  M<N 

The  number  of  nodes  can  vary  as  well  as  the  number,  length,  and 
capacity  of  the  links.  Our  problem  is  to  find  an  optimum  set  in  these 
parameters.  By  optimum  we  mean  that  design,  costing  equal  to  or  less 
than  a fixed  amount,  D,  which  maximizes  the  network  throughput  (or 
minimizes  network  GOS,  G)  subject  to  a minimum  acceptable  throughput 
after  an  attack  at  some  fixed  level. 
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2.  REPORT  ORGANIZATION 


,;__2aiS3KiH 


This  report  is  organized  in  a problem  solving  sequence.  That 
is,  we  select  a "core"  problem  from  the  previous  subsection  which 
contains  the  essential  mathematical  features  of  the  overall  problem, 
but  is  simpler  in  nature.  After  solving  the  core  problem,  we  begin 
an  augmentation  process  to  fold  into  the  core  the  remaining  features. 

Our  core  problem  is  addressed  in  Section  II.  This  problem 
involves  the  following  assumptions: 

a.  Node  costs  are  temporarily  held  constant  and  do  not  enter  into  the 
optimization  problem. 

b.  M and  N are  fixed. 

c.  No  attacks  occur  on  the  network. 

The  resulting  optimization  problem  is  then  solved.  A number  of  subissues 
occur  in  the  solution  of  this  problem  and  occupy  the  greatest  part  of  section  II. 
The  reader  may  occasionally  have  to  refer  to  the  end  of  the  section,  where 
the  subissues  are  tied  together,  to  understand  why  a particular 
subissue  is  being  addressed. 

Section  III  builds  on  Section  II  by  folding  in  the  node  cost  model 
and  allowing  M and  N to  vary.  Results  are  presented  along  the  way  to 
show  how  the  behavior  of  the  model  changes  in  response  to  the  added 
considerations. 

Section  IV  concludes  the  building  process  by  folding  in  attack 
considerations  on  the  model  of  section  III.  In  order  to  develop  these 
considerations  additional  assumptions  are  required,  which  are  given  in 
section  IV. 

Section  V discusses  the  results  of  some  typical  problems,  and 
potential  applications  and  extensions  of  the  model. 
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II.  THE  CORE  PROBLEM 


This  section  considers  the  problem  of  Section  I,  2 under  the 
simplifing  assumptions: 

a.  Node  costs  are  zero. 

b.  M and  N are  fixed. 

c.  No  attacks  occur  on  the  network. 

A number  of  subissues  are  involved  in  this  problem.  One  is  the  number 
of  tandem  links  used  in  placing  the  average  call  over  the  most  direct 
route.  This  quantity  will  be  called  L.  A second  subissue  is  the 
length,  in  miles,  of  the  average  link  in  the  network.  This  quantity 
will  be  called  -C.  These  quantities  are  somewhat  difficult  to  solve 
for  and  occupy  subsections  1 and  2 respectively. 

Subsection  3 develops  the  equations  for  a number  of  simpler 

variables  and  also  presents  a (so  far  as  we  know)  new  model  for 
computing  network  60S.  It  also  ties  all  the  pieces  together  to  solve 

the  optimization  problem  of  this  section.  Subsection  4 addresses  the 
issue  of  the  optimum  routing  strategy  for  our  type  of  network  through 
the  parameter  L. 

1 . SOLUTION  FOR  L 

L will  depend  on  a parameter  called  the  average  incidence  degree, 
I,  of  nodes  in  the  network  as  well  as  M and  N.  Let: 

I = average  number  of  links  connected  to  a node  in  the  network. 
This  parameter  tells  us  the  total  number  of  links,  T,  in  the  network. 
Since  each  link  is  connected  to  two  nodes,  and  there  are  M N nodes, 

T = MNI/2. 


To  find  L(I,M,N)  we  must  consider  particular  cases.  In  fact  we 
consider  the  cases: 

I = {2,  4.  8.  (MN-1)}.  (2) 


After  finding  these  we  assume  that  L(I,M,N)  can  be  calculated  by  Tog- 
Polynomial  interpolation  for  the  remaining  values  of  I.  Note  that  I 
cannot  be  greater  than  (MN-1)  since  at  this  value  the  network  is 
completely  connected;  i.e.,  every  node  is  directly  connected  to  every 
other  node. 

a.  Solution  for  I f2.M.N).  Figure  3 is  used  to  calculate 
L(2,M,N}.  The  problem  is  simplified  by  redrawing  Figure  3 as  a "loop" 
in  Figure  4. 

From  assumptions  A1  and  A2  we  can  look  at  any  node  in  Figure  4 and 
calculate  the  average 


L(2 


MIN{(i-l),(MN+l-^)}.  (3) 

The  MIN  (•)  operator  denotes  selection  of  the  shortest  path.  Note  that 
for  computational  purposes  we  are  assuming  that  the  node  originating 
the  call  is  labeled  1.  From  the  symnetry  of  the  problem,  (3)  can  be 
rewritten  as: 


: MN  even . 


Solving  (4)  results  in 


L(2,M,N) 


For  large  MN,  say  MN>10,  (5) 


: MN  odd 


: MN  even  . 
is  approximately 


(5) 


L(2»M.N)  = MN/4.  (6) 

b-  Solution  for  L(4.M.N).  The  geometry  of  Figure  5 applies  in  the 
solution  for  L(4,M,N).  _A  change  in  tactics  is  required  to  solve  this 
problem.  First,  we  place  an  .t,/  discrete  coordinate  system  on  the 
nodes  in  the  area.  Let 


-cl  * the  random  variable  indicating  the  -c  coordinate  of  the 
calling  node.  I = 1,2,**-,M 

yi  * the  random  variable  indicating  the  j coordinate  of  the 
calling  node.  / = 1,2,-**,N 

■c2,j2  * random  variables  indicating  the  -c  and  j coordinates  of 
the  called  node. 


From  assumptions  A1  and  A2, the  probability  density  functions  are: 

f^lU)  = = 1.2,'--,M  (7) 

fyi(j)  = fy2(y)  = VN  ; y = 1.2.*--,N.  (8) 

The  coordinate  differences  between  calling  pairs  of  nodes  are 
defined  as 

« ■ lil  - /2|  ■ (9) 
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Fitjure  5.  Geometry  for  Calculating  L(4,K,N) 


y s - X2i . 


(TO) 

Assuming  for  the  moment  that  a node  may  also  call  itself  (we  correct 
for  this  laterX  we  find  by  convolution  and  a transformation  the 
probability  density  functions 


1/N 

: y * 0 

f,(y)  = 

2(N-y)/N^ 

; y » 1,2,*'*,N-1 

(11) 

0 

: elsewhere  . 

1/M 

: ^ = 0 

f,M  = 

2(M-.t:)/M^ 

: X.  = 1,2,---,M-1 

(12) 

0 

: elsewhere 

L(4,M,N)  can  be  written  as  a function  of  the  random  variables  x 
and  y.  These  assume,  however,  that  a node  may  call  itself.  This  event 
occurs  with  probability  1/MN  and  implies  a path  length  of  zero.  To  see 
the  effect,  we  write  the  equation 


(0)  + ^ (L(4,M,N))  = E{x}  + E{y}  . 

Thus,L(4,M,N),  corrected  to  exclude  the  assumption  that  a node  may  call  itself, 
is 


L(4.M,N)  = j;^  (E{x>  + E{y}). 

The  sum  of  the  expectations  of  x and  y symbolize  the  accounting  scheme 
associated  with  finding  path  length  in  the  geometry  of  Figure  5.  Using 
(11)  and  (12)  we  find 


L(4,M,N) 


2(M-x.)  . 


(13) 


12 


i 

I 

i 


Solving  (21)  requires  use  of  the  identities  [4] 


= 

2;  ^2.n(m-l)(2KH-ll 
^»1 

2;  ^3.  slnVli^ 

^=1 

n 

= yn+l)(2n+l)(3nWl) 

^=1 

and  a lot  of  algebra  which  we  leave  to  the  ambitious  reader.  The 
solution  we  found  is 

U(8.H.N)  = 

♦ ^N2(N.l).MN(M.))-i^2M3-3M^+M)-^2N^-3N^+N)j.  (22) 

d.  Solution  for  L((MN-1).M.N).  This  case  is  simple.  Since  every 
node  is  directly  connected  to  every  other  node 

L((MN-1)^M,N)  = 1.0  . <23) 

e.  Summary  of  Solutions  for  L(I.M.N).  So  far  we  have  obtained  the 
number  of  tandem  links,  for  the  average  call,  on  the  shortest  route  for 

an  M by  N node  grid  network  for  various  node  incidence  degrees,  assuming  a 
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uniform  traffic  distribution.  We  assume  intermediate  node  incidence  degrees 
are  adequately  approximated  by  a log-polynomial  interpolation  . Several 
interpolation  schemes  were  tried.  We  finally  settled  on 

L(I.M.N)  [Log  + a2[Log  (24) 

:I  = ,MN-1  . 

The  results  so  far  are  summarized  below  and  an  example  plot  for  M=6, 

N=9  is  given  in  Figure  7. 

L(2,M,N)  * MN/4  (large  MN) 

1 f/l  M Nl  - ^ Fm^+SM-I  ^ n2+3N-1  ol 

L(4.M,N)  - pupp  ZJ 

L(8.M.N)  = - ^^(21M^-30M^+5m2+4) 

+ ^N^(N-1)-MN(M-1)  + ^2M^-3M^+m)  - 1(2N^-3N^+N)]| 

L((MN-1),M,N)  = 1.0  . 

2.  SOLUTION  FOR  I 

Our  core  problem  assumes  that  only  links  in  the  network  cost 
money.  To  optimize  the  design  we  must  minimize  the  network  GOS  by 
distributing  the  channel  miles  of  links  we  can  purchase  in  an  optimum 
way.  The  first  thing  we  must  discover  is  how  many  link^  we  can  buy; 
this  is  T in  equation  (1).  Next  we  need  to  know  the  capacity,  c,  of 
the  average  link.  To  find  this  we  must  know  since  the  amount  to  be 
spent,  0,  is,  from  A4: 

D = T feci  . (29) 

+ This  report  used  10  as  the  log  base  throughout.  Any  other  base  will  work 
equally  well  as  long  as  consistency  is  exercised. 
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(25) 

(26) 

(27) 

(28) 


Jjl 


Thus  1f  we  know  we  can  compute  c. 


As  in  the  previous  subsection,  we  must  consider  particular  cases  and 
the  dependence  of  ^ on  I (as  well  as  M and  N).  The  cases  considered 
are 

I = {2,4,8.(MN-1)}  . (30) 

After  finding  these  we  assume  that  £(I,M,N)  can  be  calculated  by  log- 
polynomial  interpolation  for  the  remaining  values  of  I. 

One  of  the  factors  involved  in  computing  t is  the  distance,  in 
miles,  between  nodes  at  (-c,/)  and  (-c+l,/);  or  equivalently  (-c,y)  and 
U»y+1).  This  is  our  unit  distance  u.  Since  we  are  treating  a square 
grid  type  of  network  with  MN  nodes  covering  an  AB  square  mile  area: 

a = [AB/MN]^^^  . (31) 

a-  Solutions  for  •e(2.M,N) , -£(4,M.N),  and  -^(8.W,N.  The  following 
equations  can  be  written  by  inspection  from  Figures  3,  5,  and  6 and 
equation  (31) 


£(2,M,N)  = [AB/MN]^/^ 

(32) 

£(4,M,N)  = [AB/MN]^^^ 

(33) 

^(8,M,N)  = [AB/MN]^/^  . 

(34) 

b.  Solution  for  £((MN-1),MN).  The  solution  for  -eClMN-U.M.N)  is 


not  quite  obvious.  This  situation  corresponds  to  a completely  connected 
network,  however,  and  some  thought  reveals  that  it  is  mathematically 
analogous  to  the  problem  of  finding  the  straight-line  distance  between 
calling  pairs  of  nodes  for  the  average  call. 


As  a consequence  of  this  analogy,  we  can  write  £((MN-1),M,N)  in 

MN 

terms  of  (11),  (12)  and  (31),  and  the  correction  factor  jjjj^  first 
used  in  (13).  That  is. 


£((MN-1),M,N)  = (AB/MN)^/^  E{(x^^y^)^/^} 

or 


^((MN-1),M,N) 


(35) 


M-1 

I- 

-1=1 


m2 


M-1 

-i=l 


N-',, 2 2)1/2  4(H-^)(jha 

i=l 


(36) 


This  reduces  to 


£((MN-1),M.N)  = (AB/MN)^/^ 

+ X X (^V)^^^  4(M-^j(N-jl1 

-1=1  /=!  J ■ 


(37) 


c.  Summary  of  Solutions  for  £(I,M,N).  In  this  subsection 
(11,2),  we  have  obtained  the  average  link  length  for  an  M by  N 


node  grid  network  for  various  incidence  degrees.  We  assume  intermediate 
incidence  degrees  are  adequately  approximated  by  a log-polynomial  inter-, 
polation.  After  trying  several  interpolation  schemes,  we  finally  settled  on 

^(I,M.N)  =^0+0^  [Log  |]  + a2[Log  (38) 

: I = 4,5,  ...MN-1. 


The  results  are  summarized  below,  and  an  example  plot  for  M=6,  N=9  is  given 
in  Figure  8. 


^(2,M.N)  = [AB/MN]^/2 
^(4,M,N)  = [AB/MN]^/^ 

^(8,M,N)  = [AB/MN]^/^  [■^“^1 
£(MN-1),M,N)  = [AB/MN]^/2[^^m2+n2-2) 

- 2 2uV)’/2 


(39) 

(40) 

(41) 


(42) 


3.  PUTTING  IT  TOGETHER 

The  reader  has  undoubtedly  noticed  by  now  that  a key  parameter 
in  the  subissues  of  our  core  problem  is  the  node  incidence  degree,  I. 
This  observation  suggests  an  algorithm  which  computes  network  GOS  as 
a function  of  I.  Suppose,  then,  that  we  are  given: 

A : width,  in  miles,  of  the  rectangular  area  to  be  serviced 
B : length,  in  miles,  of  the  rectangular  area  to  be  serviced 
D ; the  amount  to  be  spent  on  the  network 
E : network  offered  traffic  in  erlangs 


M : number  of  nodes  along  the  A dimension  of  the  rectangular  area 
N : number  of  nodes  along  the  B dimension  of  the  rectangular  area 
I : the  average  node  Incidence  degree 
fe  : cost  per  channel -mile  for  links. 

Now  we  can  compute,  from  (1) 

T = MNI/2  (43) 

: total  number  of  links  In  the  network. 

Next  we  can  find,  from  (43),  (29),  and  (38)  through  (42) 

c = D/Tfce 

= 2D/MNIfe£  (44) 

: the  capacity  of  the  average  link. 

Note  that  I = £(I,M,N)  Is  found  by  the  procedure  of  section  11,2-3. 
By  our  assumptions  of  uniform  traffic  distribution  and  network  synrietry 
we  can  focus  our  attention  on  a "typical"  node  in  the  network.  This 
node  1s  characterized  ty  Its  average  incidence  degree,  I,  the  average 
link  capacity,  c,  and  the  originating  traffic  at  that  node,  Eq;  i.e., 

= E/MN  . (45) 

Suppose,  for  the  moment,  that  network  GOS,  G,  is  some  function  of 
c,I,Eq,  and  L(I,M,N);  i.e., 

G = g(c,I,Ep,L)  (46) 

where  L = L(I,M,N)  is  given  by  the  procedure  of  section  II,  l,e.  Then 
our  core  problem  is  solved  by  the  algorithm  of  Figure  9. 


a.  Solving  for  G.  By  network  GOS,  6,  we  mean  the  probability  that 
the  average  call  is  blocked  by  the  network.  Thus  an  equivalent  measure  is 
network  call  completion  rate,  R,  which  is  the  probability  that  the  average 
call  is  completed: 

R = 1.-6  . (47) 

A number  of  techniques  are  available  for  calculating  these 
measures,  in  an  approximate  way,  for  the  more  traditional  discrete 
network  design  problem.  None  are  directly  suitable  for  our  purposes, 
for  two  reasons: 

a.  All  depend  on  an  assumption  of  statistical  independence  of 
traffic  from  one  link  to  any  other  - which  is  certainly  not  a 
valid  assumption. 

b.  All  take  a "microscopic"  view  of  the  discrete  network  structure; 
whereas  our  approach  is  "macroscopic"  in  nature  and  deals  with 
average  properties  of  the  network. 

Consequently  we  develop  an  approach  which  is  felt  to  be  similar  in 
spirit  to  the  Katz  algorithm  [5],  but  calculates  R from  first  principles 
for  our  core  problem. 

To  begin  with  we  need  to  observe  that  a node  operates  on  three 
identifiable  categories  of  traffic: 

Eq  ; originating  traffic,  in  erlangs. 

This  is  calling  traffic  originating  in  the  access  area 
services  by  the  node. 
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Ep  : destination  traffic,  In  erlangs. 

This  Is  traffic  which  originated  elsewhere  In  the  network  and 
Is  destined  for  users  In  the  access  area  serviced  by  the  node. 

Ej  : tandem  traffic.  In  erlangs. 

This  Is  traffic  which  originated  elsewhere  In  the  network  and 
Is  destined  for  users  not  In  the  access  area  serviced  by  the 
node. 

From  these  definitions  and  assumptions  A1  and  A2  It  follows  that 
R - Ep/E^  . (48) 

We  need  to  find  Ep. 

Ep  must  comprise  traffic  entering  our  "typical"  node  over  links 
incident  on  that  node.  However,  not  all  of  the  traffic  entering  a node 
belongs  to  Ep.  In  fact  we  need  a new  parameter  to  describe  a2i  traffic 
entering  a node  via  backbone  links.  Call  this  parameter  Ej.  By  the 
symmetry  (average  properties)  of  the  network  in  question,  as  much  traffic 
is  flowing  out  of  the  node  onto  the  links  as  is  flowing  into  the  node  from 
the  links. 

Ej  : the  sum  of  all  traffic,  in  erlangs,  entering  a node  from  the  . 
links  incident  on  the  node.  Conversely,  It  Is  also  the  sum  of  all 
traffic.  In  erlangs,  leaving  a node  and  entering  links 
Incident  on  the  node. 

Again,  by  symmetry,  we  see  that  this  parcel  of  traffic,  Ej  is  divided 
equally  amongst  all  links  Incident  on  the  node.  Thus  each  link 
carries  E|^p  erlangs  where 

Elc  = E^'/I  . (49) 

26 


How  much  traffic  Is  offered  to  each  link?  Call  this  quantity  E^^q.  Then, 
by  symmetry 

Elo  = (Eq+Ej)/!  • (50) 

Now,  to  a first  approximation,  Ej^q  and  Ej^j,  are  related  by  the  Erlang  B 
equation 

^LC  ' (^^^ 

q = Eg(E^Q,c)  (52) 

and  Eg(‘,*)  denotes  the  Erlang  B equation. 

We  need  one  more  observation;  this  concerns  the  origin  of  Eg. 

Eg,  by  the  definitions,  must  be  some  fractional  part  of  Ej.  The  remaining 
fractional  part  of  Ej  must  be  Ej.  The  preceding  discussion  is 
summarized  in  the  block  diagram  of  Figure  10. 

We  shall  write  an  equation  for  solving  Figure  10  presently;  but 
first  we  need  to  determine  the  fractional  relationship  between  Ey  and  Eg. 
At  this  point  L(I,M,N)  becomes  useful.  Suppose,  for  the  moment 
that  L(I,M,N)  represented  the  true  average  number  of  tandem  links 
per  call,  rather  than  that  average  over  the  shortest  route.  Now  if  we 
randomly  select  a traffic  parcel  on  a link, we  need  to  determine  the 
probability  that  that  parcel  exits  the  network  (becomes  part  of  Eg)  at 
the  node  being  fed  by  the  link.  This  line  of  thinking  shows  that 

Eg  = Ef/L(I,M,N)  . (53) 


From  (48)  and  (53)  we  have 
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Note  that  (62)  is  transcendental  due  to  folding  in  (61).  We  can  solve 
(62)  for  Ej,  and  thereby  find  R and  G via  (55)  and  (47)  by  using  a 
simple  search  algorithm  which  we  need  not  go  into  here.  Suffice  it  to 
say  that  this  search  algorithm  performs  the  function 

G(I)  = 3(c,I.Eq,L(I.M,N))  (63) 

in  Figure  9.  One  fine  point  needs  to  be  mentioned.  That  is,  due  to 
the  averaging  properties  used  herein,  c is  not  necessarily  integer  in 
(61).  Therefore,  we  use  a logarithmic  interpolation  technique  suggested 
in  [7],  for  finding  Eg  (p,c)  when  c is  not  integer. 

b.  Some  First  Results  on  the  Core  Problem.  We  must  still  address 
our  earlier  assumption  that  L(I,M,N)  represents  the  true  average  number  of 
tandem  links  per  call.  This  is  done  in  the  next  subsection.  First,  however, 
we  use  that  assumption  to  show  the  results  of  a problem  solved  with  the 
algorithm  of  Figure  9. 

The  problem  parameters  are: 

M=6  N=9 

A = 2000  miles,  B = 3000  miles 

D = $2,000,000  E = 4800  erlangs 

fe  = $0.50  per  channel  mile. 

The  results  are  shown  in  Figure  11. 

With  this  graph  one  can  pick  out  I for  G*,  the  minimum  GOS,  and 
from  this  value  of  I,  I*,  L(I*,M,N),  and  £(I*,M,N)  can  be  determined. 

The  average  link  capacity,  c,  can  be  computed  through  (1),  (28),  and 
Figure  8. 


4.  R DEPENDENCE  ON  L(I.M,N) 

Wp  promised  earlier  to  investigate  the  assumption  that  L(I,M,N) 
represents  the  true  average  number  of  tandem  links  per  call.  It  clearly 
is  the  minimum.  But  what  if  alternate  routes,  using  more  tandem  links 
than  the  shortest  route,  are  allowed  for  calls?  The  net  result  will  be 
to  increase  the  average  nutrtber  of  tandem  links  per  call.  Some  people  argue 
that  this  effect  increases  network  GOS  by  increasing  the  effective  loading 
on  the  network;  others  will  claim  that  network  GOS  is  reduced  since  each 
call  now  has  additional  chances  of  being  placed  in  the  network.  Figure 
10  gives  us  an  opportunity  to  mathematically  investigate  this  issue  for 
our  type  of  network. 

All  the  arguments  of  section  II,  3, a hold  now  as  before.  The 
difference  is  that  L(I,M,N)  must  be  replaced  by  L'(I,M,N)  and 

L'(I,M,N)  > L(I,M,N)  (64) 

Since  L(I,M,N)  represents  the  minimum  value  on  tanderp  links  for  the 
average  call. 

What  we  must  do  then  is  investigate  R in  (55)  for  its  behavior  as 
L'  increases.  We  drop  the  (I,M,N)  argument  on  L and  L'  for  the  remainder 
of  this  subsection.  From  (55) 

R = MN/EL'  = f(EpL'  ) (65) 

and  we  must  find 
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r 


t 

I 


^ t L'  ).  (66) 

At  this  point  we  can  observe  that  if  ^ is  everwhere  negative  for  L > L» 
it  Is  established  that  the  model  of  section  II,  3, a is  valid  as  is;  and 
we  can  conclude  that  use  of  shortest  path(s)  only  produces  optimum 
network  performance.  If  this  is  not  the  case,  then  we  can  find  an  L*, 
which  is  not  the  same  as  L,  and  which  produces  a maximum  network  call 
completion  rate. 

Equation  (66)  can  be  attacked  by  using  the  chain  rule  for 
differentiation  of  composite  functions  [8] 


dR  _ 3f(E;,L)  dE;  8f(E'  L') 


(67) 


Now 


If  =[f  - ^ Eg 6t)  c') 
,|[i  - Ee(Ml  -,xOil}c)] 


and 


*('t)  ['-Eg  (MlirillM.  • 


If  we  make  the  following  identifications 
'E; 


-n  E'  0 
“ Ll'  ■ ■ mJ 


^1  ■ ii'h 


S = E ( • *1 

^2  aEf  ""B  ^ ^ 

P = 1 - Eg  (•.  •) 


Then  (71)  becomes 


Ey 


dE^  ^^1  l'2  ^ 


xsg  + (i4)p 


Using  the  relation  [4] 


d r \ T 


(72) 


(73) 

(74) 

(75) 

(76) 

(77) 


(78) 


(79) 
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and  letting 


"[wi  Ot'3 

we  see  that 

and 

^2  ^B  (’"*c)  j- (l-p/ 

Now  if  we  Identify 
^3  ' ^ 


then  (78)  becomes 


Combining  this  with  (68)  and  (69)  Into  (67) 


Et  MN 


dR  ^ L'  \ ^ , 

Clearly 


^ < 0 for  > L>1  . 


(85) 


(86) 


This  result  proves  two  things  for  our  network  problem: 

a.  Maximum  network  call  completion  rate,  R,  (or  minimum  network 
60S,  G)  is  obtained  when  shortest  path(s)  only  are  allowed  in 
placing  calls. 

b.  From  this  it  follows  that  the  model  of  Section  II,  3, a for 
calculating  R is  valid,  for  computing  optimum  R,. as  Is.. 


III.  SWITCHES  COST  MONEY  TOO 


The  purpose  of  this  section  is  to  augment  the  core  problem  of 
section  II  by  folding  in  the  node  cost  model  given  by  assumption  A5. 

1 . FOLDING  IN  THE  NODE  COST  MODEL 

By  assumption.AS, the  cost  for  a single  node  is 

Dj,  = a + fat  . (87) 

Now  assume  for  the  moment  that  M,  N,  and  I are  fixed.  The  link  cost 
model  given  by  A4  is 


= kcZ  . (88) 

Since  M,  N,  and  I are  fixed,  the  total  number  of  links  is 

T = MNI/2  . (89) 

From  these  we  see  that  the  amount  spent  on  all  nodes,  Dj|^,  is 

= [a  + faIc]MN  (90) 


and  the  amount  spent  on  all  links,  D-.  , is 

I 


° °TL 


(92) 


Substituting  (90)  and  (91)  Into  (92)  and  solving  for  c 

c » [D-aMN]/[6MNI  + feMNI^2]  (93) 

where  t = £(I,M,N).  Using  (93)  instead  of  (44)  in  the  algorithm  of 
Figure  9 accounts  for  node  costs. 

2.  ALLOWING  M AND  N TO  VARY 

Allowing  M and  N to  vary  in  the  modified  (by  (93))  algorithm  of 
Figure  9 requires  only  the  addition  of  an  outer  loop  to  step  through 
various  values  for  M and  N.  Two  observations  are  worth  making: 

a.  We  required  early  on  that  the  nodes  be  arranged  in  a square 
grid.  In  the  spirit  of  this  it  is  necessary  to  insure  that 

I • B • (M) 

b.  Changing  the  number  of  nodes  in  the  complete  A x B area  does  more 
than  change  the  backbone  network  behavior.  It  also  impacts:  (1) 
the  number  of  access  areas  (one  per  node),  and  (2)  the  cost  of 
"wiring"  the  users  in  the  access  area  to  the  node  connecting  them 
into  the  backbone.  For  this  report  we  ignore  the  second  issue. 

It  can  (and  will)  be  incorporated  into  our  model  at  a later  date. 

3.  SOME  RESULTS 

Figure  12  shows  the  results  of  a run  using  the  augmented  algorithm 
of  Figure  9.  Two  observations  can  be  made  from  this  run: 

a.  Network  performance  is  steadily  improved  as  the  number  of  nodes 
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Is  reduced.  This  effect  is  in  agreement  with  intuitive 
"econony  of  scale"  arguments  since  it  implies  links  with 
larger  capacity. 

b.  The  node  incidence  degree  for  optimum  network  performance  is  not 
as  sharply  defined  here  as  in  the  earlier  case  where  nodes 
were  "free".  Compare  Figure  12  with  Figure  11.  The  effect 
here  seems  to  be  explained  by  costs  of  terminating  channels 
onto  a node.  Fewer  terminations  are  required  if  the  incidence 
degree  increases,  since  this  causes  average  link  length  to 
increase.  In  turn, average  tandem  links  per  call  decrease, 
which  tends  to  improve  network  60S. 
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IV.  SURVIVABLE  TOPOLOGIES 


The  preceding  sections  have  dealt  with  the  problem  of  minimizing 
network  GOS  under  constrained  costs.  It  was  assumed  that  no  attacks 
occurred  on  the  network.  Thus  the  network(s)  found  were  efficient 
under  benign  conditions.  When  designing  the  DCS  it  is  also  necessary 
to  characterize  the  networks s)  under  attack  conditions. 

1.  THE  PERFORMANCE  CHARACTERISTIC 

Reference  [1]  develops  the  idea  of  "Performance  Characteristic"  in 
detail.  Here  we  present  only  the  essential  features.  Let  a family  of 
network  designs  be  characterized  by  (M,N).  Now  let  the  most  efficient 
design  in  that  family,  under  benign  conditions,  be  characterized  by 
(I*.M,N);  where  I*  is  that  node  incidence  degree  which  minimizes  network  GOS. 
The  algorithm  of  Figure  9 can  be  used  to  find  (I*,M,N)  and  produce 
G(I*,M,N)  for  that  design.  We  define  the  utility  under  benign 
conditions,  Uj^,  as 

Ujj  = E(1-G(I*,M,N))  . (95) 

With  that  same  design  (I*,M,N),  we  associate  a quantity  called  the 
utility  after  an  attack,  U..  U,  is  found  by  removing  from  that  design 

a a 

a specified  number  of  nodes,  w,  and  calculating  the  "throughput"  of 
the  resulting  network.  That  is 

Ua  = (1-G(I*',m',n'))  (96) 


41 


where 


= MN-W  . (97) 

The  following  assumptions  are  implicit  in  (96); 

A8.  When  a node  is  removed  from  the  network,  the  traffic  in 
the  access  area  serviced  by  that  node  is  lost. 

A9.  No  cuts  occur  in  the  network. 

AlO.  The  post-attack  network  still  has,  on  the  average,  a square 
grid  structure  characterized  by  M,'  N,'  and  I*C  where  H'  N^is  given  by 
(97)  and  I* 'Is  developed  as  shown  below. 

The  pre-attack  network  has  a total  number  of  links,  T,  given  by 

T = MNI*/2  . (98) 

Since  W nodes  are  removed  to  develop  the  post-attack  networks,  and 
each  node  has,  on  the  average,  I*  links  incident  on  it;  the  total  number 
of  usable  links,  T^  in  the  post-attack  network  is 

t'  = (MNI*/2)  - WI*  . (99) 

There  are  (MN-W)  = M^N^  nodes  in  the  post-attack  network.  Thus  its 

( 

average  incidence  degree,  I*  , is 

l/  = (MNI*-2WI*)/(MN-W)  . (100) 

The  implicit  assumption  is 

All.  Nodes  removed  in  the  attack  are  not  adjacent  to  each  other. 
The  preceding  conditions  allows  use  of  the  model  of  section  II,  3, a for 
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I finding  The  necessary  relationships  are 

a.  Redefining  the  post-attack  network  offered  load,  E , as 

I 

E = E(MN-W)/MN  (101) 

In  accordance  with  assumption  A8. 

b.  Computing,  for  (M^,N^),  the  array  L(I,M',N'’)  and  finding  from 

this  the  value  of  L(I*',M^n').  ^ 

The  results  of  this  subsection  characterize  the  most  efficient  design  ; | 

in  a family,  (I*,M,N)  by  a pair  of  numbers,  and  Ug.  Clearly  we  I 

can  now  repeat  this  process  by  varying  M and  N to  trace  out  a locus  of  j 

points  (Ujj,Ug)  in  two-dimensional  space.  This  locus  of  points  is  the 
performance  characteristic. 

2.  SOME  MORE  RESULTS 

The  preceding  subsection  developed  a new  way  of  thinking  about 
how  well  a network  design  performs,  j Some  results  are  given  in 
Figures  13  and  14.  As  the  figures  show,  the  presentation  is  similar 
to  the  "guns  versus  butter"  curves  seen  in  economics  textbooks.  For 
^ example.  Figure  13  shows  that  we  can  achieve  almost  non-blocking 

performance  under  benign  conditions  by  building  a backbone  terrestrial 
network  with  4 x 6 nodes.  Unfortunately,  the  performance  of  that 
network  after  an  attack  which  removes  nine  nodes  is  very  poor. 

The  presentation  format  allows  an  additional  consideration  - 
, survivability  - to  be  folded  into  network  design  from  the  first  stages 

• of  design.  That  is,  we  can  specify  in  advance  a minimum  acceptable 

post-attack  throughput  for  the  contemplated  design.  Thus,  if  the 

I ! 

13 


LEGEND:  (•,«,•)  = (I*,M,N) 


DATA: 

Figure  1 


A=2000  B=3000  0=4500000  E=4800 
a=8000  b=106  k=0.5 

3.  Performance  Characteristic  for  a Grid  Network;  W=9 


LEGEND;  (•.•,•)  = (I*,M,N) 


DATA:  A=2000  B=3000  0=4500000  E=480D 

a=8000  b=106  k=0.5 

Figure  14.  Performance  Characteristic  for  a Grid  Network;  l!=10 


requirement  was,  say,  3000  erlangs  throughput  after  an  attack  removing 
nine  nodes,  the  network  design  must  have  on  the  order  of  7 x 9 nodes  (as 
shown  in  Figure  13)  with  an  average  incidence  degree  of  about  12. 


V.  CONCLUSIONS 


Approximation  techniques  developed  in  this  report  permit  analyzing 
average  structural  properties  of  terrestrial  backbone  networks.  In 
addition,  a new  way  of  measuring  the  performance  of  a network  design 
allows  survivability-capacity  tradeoffs  to  be  made  early  in  the  design 
process  where  such  decisions  are  most  crucial. 

The  concepts  introduced  in  this  report  can  be  developed  in  more 
detail  for  implementation  in  traditional  discrete  network  design 
algorithms  if  desired.  However,  the  analysis  of  average  structural 
properties  may  be  sufficient  to  guide  future  DCS  planning.  In  this  case 
serveral  extensions  to  this  report  are  desirable: 

a.  Access  area  costs  need  to  be  considered  in  the  model. 

b.  Satellite/terrestrial  tradeoffs  must  be  considered. 

c.  Integration  issues  must  be  addressed. 

d.  Various  switching  possibilities  (circuit,  packet,  SENET  [9]), 
must  be  analyzed. 

These  issues  represent  not  only  a prospectus  for  future  work,  they  also 
indicate  the  range  of  possible  applications  of  models  using  the  techniques 
developed  herein. 
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